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Abstract 

We apply the uniform asymptotics method proposed by Bassom, Clarkson, 
Law and McLeod [1] to a special Painleve V equation, and we provide a simpler 
and more rigorous proof of the connection formulas for a special solution of the 
equation, which have been established earlier by McCoy and Tang via the isomon- 
odromy and WKB methods. 
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1 Introduction 


We apply and extend the method of nniform asymptotics proposed by Bassom, Clarkson, 
Law and McLeod in [1] to a special case of the hfth Painleve eqnation (PV) 
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which appears in the time-dependent correlation fnnctions of the transverse Ising chain at 
the critical value of the magnetic held Oils]. Our main focus will be on the connection 
formulas of the equation. 

The equation fll.ip is reducible to a special case of the third Painleve equation (PHI). 
Indeed, if we put y{x) = with x = At, then w{t) satishes Painleve III, namely 
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which arose in Ising model studies [16]. McCoy, Tracy and Wu [16] derived oo -H- 0 
connection formulas for a two-parameter class of bounded solutions w(t;Q,p) of the 
one-parameter family of Panleve III equations fll.2p . 

The special Painleve V equation flTlh . or the equivalent Painleve III equation fll.2p . 
plays a crucial role in problems related to random matrices and random processes, or¬ 
thogonal polynomials, string theory, and in exactly solvable statistical mechanics and 
quantum field models. For example, for the special case 0 —)■ 0, if we set w{t) = —e*"^*-*^ 
in (11.21) . then 0(t) satisfies the following equation 

— 2tsin20(t) -|- 2sin0(t) = 0, (1.3) 


which appeared in the problem of classification for rotation surfaces with harmonic in¬ 
verse mean curvature [3]. This equation has been numerically investigated by Bobenko 
et al. [3], and also studied by Andreev and Kitaev [1] based on the results obtained in 
[2] by using isomonodromy deformation and the WKB method. 

The problem concerning one particle density matrix of impenetrable bosons at zero 
temperature (see Creamer et al. [5], and Vaidya and Tracy [25]) has been reduced in [5] 
to a study of the equation 

cp"{z) = (((f'(z)f - 1) cot (f(z) + -(1 - (f'(z)). (1.4) 

'Z/ 

Creamer et al. [5] have studied the solution of the equation numerically. An analytical 
study of (11.4p has carried out by Suleimanov [23] . 

It is worth pointing out that the equation (11.41) is equivalent to equation (11.3p via 
the transformation of the dependent variables e**^*^^) = i tan ^ and the independent 
variables z = 2it. 

In another special case when 0 = |, we may take the transformation w{t) = e*’^*'*^ 
and then the Painleve III equation (II.2p becomes 

H —= 2sin2'0(t). (1.5) 

iy 

For a numerical study of the last equation, see Lamb with a connection to the tt- 
pulse of the sine-Gordon equation. Flaschka and Newell [8] have also considered equation 
(II.5p by using the methods of monodromy preserving deformation and singular integral 
equations, and obtained representation of the one-parameter family of solutions to (11.51) 
that are holomorphic at the origin. 

The same differential equation (II.5p . adapting the new variable s = 4t^, is a special 
case of the following equation 


r+ 

s 
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sin(2V') 


cos-^ 
4s2 sin^ '0 


with d = 0, which is closely related to the asymptotics of Bessel kernel limit of the 
Fredholm determinant describing the statistics of the level spacing of the eigenvalues of 
Hermitian matrices of large order in a single interval; see Tracy and Widom [24] . 
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Now we see that the Painleve V fll.ll) is related to various nonlinear equations with 
statistic physical backgrounds. Yet the main objective of the present paper, is to calculate 
the asymptotics of a special solution as as x —?■ 0 and x — )■ +oo, and to justify connection 
formulas between parameters involved in the asymptotic approximations. It is known 
that there exist solutions to (11.Ih . regular on the positive real line, with behavior at the 
origin 

y{x) = ^ [l + (1 — 20)x + O {x"^)] as a; —)■ 0; (1.6) 

cf. [T51[T7] . where the parameter p is not on the imaginary axis, and such that |Im p\ ^ 0. 

McCoy and Tang m [T8] derived +oo -H- 0 and ±ioo -H- 0 connection formulas 
for two-parameter solutions of the one-parameter family of Painleve V equations fll.ip . 
respectively. In this paper, we focus on the -|-cxd -h- 0 connection problem. It is known 
from that the solution of fll.ip satisfying boundary condition (II.6p possesses the 

following asymptotic expansion 

y{x) = —1 + 4:X~^Fi{s) + 4:X~^ [20 — 1 — 2F‘^{s)] + O as x —)> -|-oo, (1.7) 

where 


Fi(s) = ae -f be 


^—2is 


( 1 . 8 ) 


with s = | —a&ln |. Here a and b are independent of x, and |Im{a6}| < j. The following 
results state that the parameters a and b in (11.81) are explicit functions in the parameter 
p in ([T6P : 


Theorem 1. 


a{p) = 7 ^P (1 + 22 c(p))e^e-®*"(^)'" 2 e""(^) 

2y TT 


^(P) = y^r(l-2ic(p))e 4h 

2 -y TT 


-Y«>6*c(p) In2g7rc(p) 


COS 7r0 — ^ sin 7r0 
cos 7r0 -|- ^ sin 7r0 


where 


(1.9) 

( 1 . 10 ) 




1 + ^1 sin^ t:Q 


( 1 . 11 ) 


Remark 1. From fll.Qp and fll.lOp it is readily observed that 

a{.p)b{p) = c(p). (1.12) 

Here use has been made of the fact that P(1 -(- 2ic)P(l — 2ic) = From fll.lip 

and fll.l2p . we see that |Im{a6}| < | if and only if p is not on the parts of imaginary 
axis with |Imp| ^ 0. 
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The above connection formulas have been established by McCoy and Tang in ini by 
using the method of isomonodromy deformation and the matching of WKB solutions. 
To explain their approach, we briefly outline the relation of Painleve V with the theory 
of monodromy preserving deformations of linear ordinary differential equations with 
rational coefficients. The reader is referred to mm for more details. 

The Lax pair of the fifth Painleve equation, with parameters ©o = 0i = 0 and 
©oo = 0, is a system of linear differential equations for the matrix function \h(A,x), 
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(1,13) 


(1.14) 


where y{x), v(x) and u(x) satisfy the following isomonodromy deformation system: 



xy - 2v(y - if -h 2e(y - 1), 


dn 2 
X— = yv 
dx 


i(t> + e)^ 
y 


x—hm = —2v — Q + yv ^—(n + 0). 
dx y 


(1.15) 

(1.16) 
(1,17) 


Furthermore, y{x) solves the fifth Painleve equation fll.ip . In this sense, the Painleve V 
equation is equivalent to the compatibility condition = fl/xA- 

The canonical solutions = Tfc(A) of fll.131) . A; = 1,2, are determined by their 
asymptotic approximations. In a neighborhood of the irregular singular point A = oo, 
they have the following asymptotic expansion: 



|A| —)■ oo, X e M"*", /c = 1,2, 


respectively in the Stokes sectors fli and 112 defined as 


(1,18) 


Hfc 



+ ■n{k — 2) < arg A < 


Stt 

Y 


+ 7r(A: — 2), |A| > i? 


A: = 1,2, 


with arbitrary finite positive constant R. 

These functions are related by certain Stokes matrices Si and 5*2, namely, 

(1.19) 

where 
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and the constants si and S 2 are called the Stokes mnltipliers. 

The canonical solutions at A = 0 and A = 1 are 2x2 uniniodular matrices \1/^‘’^(A, x) 
and \1/(^)(A, x), and can be dehned as 

^5°°)(A,a:) = T(°)(A,x)Eo, ^S°°Ha,x) = a:)Ei, (1.20) 

where Eq and Ei are unimodular constant matrices. Furthermore, the behavior at the 
singularities are 

¥^\X,x) = asA^O, (1.21) 

T(^)(A,x) = T(^)(A,x)(A- 1)^®“"^ asA^l, (1.22) 

with T®(A,a:) and T*'^^(A,a:) being holomorphic at A = 0 and A = 1, respectively. 

The monodromy matrices E^ at the regular singularities X = k {k = 0,1) fulhl the 
following cyclic relation: 


E-i^-^ieasE^ = 


(1.23) 


which specihes Eq and Ei up to left-multiplicative diagonal matrices diag((io, <^o ^) and 
diag((ii, d¥), respectively. The connection matrices Eq and Ei are independent of x\ The 
isomonodromic condition = 0 holds here owing to the fact that all of \1/^°°^(A, x), 
'l/(°^(A,x) and T*^^)(A,x) are solutions of fll.ldp . We note that none of the entries in 
matrices Eq and Ei vanishes. In fact, for example, if (i?o) 2 i = 0, then from [2] we see 
that all the Stokes multipliers = 0, fc = 1, 2. From the results in HDl for the connection 
matrices and in view of fll.23p . we obtain p = iQ, which contradicts the restriction on p. 

Appealing to the isomonodromic deformation techniques developed in [9], we may 
use Eq and Ei to express the parameters a and b in (II.7p as functions of the initial 
parameter p in fll.bp . More exactly, the quantities to be used are 


{Eq)ii{Eq)22 

{Eq)i2{Eq)21 


and 


(A'i)ii(A'i)22 

(-E'i)i2(-E'i)21 


(1.24) 


(1.25) 


As mentioned earlier, the connection problem fll.9p - fll.10p was hrst established by 
McCoy and Tang d- They considered the asymptotic behavior of the solution \k(A,x) 
to the hrst order system fll.l3p as a: —)■ oo for A being kept away from the turning 
point i and the singular points 0,1. Then, they matched these WKB solutions with the 
asymptotic approximation at the turning point involving parabolic cylinder functions, 
and the behavior at the singularities A = 0,1 involving the Bessel functions. Eventually 
they were capable of calculating the two invariants fll.24p and fll.251) as x —)■ -|-oo. The 
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invariants Jq and Ii for small x can be obtained from the connection matrices Eq and 
El, which have been done by Jimbo [10]; see also diiiH], such that 


lo 

h 


ip + Q 
ip-e' 
ip — Q 
ip + Q' 


and 


(1.26) 

(1.27) 


The facts that Jq and Ji are independent of x will lead to equalities involving a, b and 
the initial parameters. The formulas in Theorem [T] would then follow. The procedure, 
however, is complicated, and is difficult to make rigorous; see a comment made in jH 
p.245]. 

In this paper, we shall provide a hopefully simpler and more rigorous derivation of the 
formulas fll.Op - fll.lOp by using the method of uniform asymptotics proposed by Bassom, 
Clarkson, Law and McLeod |1]. Along the same lines we may hnd the work of Olver 
[20] and Dunster [B] for coalescing turning points. Initially in |1], PII has been taken as 
an example to illustrate the method. While the difficulty in extending the techniques 
for PII to other transcendents is also acknowledged by the authors of [U p.244]. Yet the 
method has been applied to the connection problems for PHI (Sine-Gordon) and PIV; 
cf Wong and Zhang [23 |28], and has also been used to hnd the asymptotic behavior at 
inhnity of the solutions to PIV [22] and PV [El [13] . 

We briehy outline the uniform asymptotics approach to derive the connection formu¬ 
las fll.9p - fll.10p . First, we will obtain the second-order differential equation fl2.6p from 
the Lax pair fll.l3p . Then, we substitute the known large-x asymptotic behaviors of 
y{x) and v{x) into the second-order equation and obtain an approximate equation (12.9p . 
The equation fl2.9p has only two coalescing turning points which coalesce with a reg¬ 
ular point. Thus, uniform asymptotic solutions are to be constructed in terms of the 
parabolic cylinder functions, uniformly for A on the Stokes curves as x —)■ cxo; cf. Olver 
[20] . and see Theorem [2] in Section [2] below. 

Accordingly, we will calculate the connection matrix via fll.20p by using the asymp¬ 
totics of the fundamental solutions on the Stokes curves; cf. Theorem [21 respectively as 
|A| —oo and A —?■ 0 (or A —)• 1), for x —)■ -|-oo. Consequently, we can use fll.24p - fll.25p to 
obtain Jq and A for large x] see fl3.52p - fl3.53p below. Finally, equating fll.26p and fl3.52p 
gives b = b{p), and equating fll.27p and fl3.53p gives a = a{p). 

The difference between the method of uniform asymptotics and the WKB method 
is that the latter has to match different approximations in different regions while in 
the uniform asymptotics cases the complicated matching procedure is not needed. The 
derivation of the method of uniform asymptotics is rigorous and may lead to a simpler 
argument with minimal computational efforts. However, for the uniform asymptotics 
method, difficulties may arise in describing the geometry of the Stokes curves near the 
turning points in the hnite plane, such as determining of the correspondence between 
±CK and Ai ,2 in fl2.15p . 

The rest of the paper is organized as follows. In Section [21 we derive uniform ap- 
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proximations to the solutions of the second-order differential equation obtained from the 
Lax pair fll.lSp as x —?> -|-oo by virtue of the parabolic cylinder functions on the Stokes 
curves. The last section is devoted to the evaluation of the two invariants Jq and Ji in 
fll.24p - fll.25h for large x. The proof of Theorem [T] is also provided in that section. 


2 Uniform asymptotics as x ^ +(X) 

In the present section, we apply the uniform asymptotic method to deal with the large- 
X behavior of the second-order differential equation fl2.9h obtained from the Lax pair 
([ITS]). 

First, we eliminate the function u{x) in fll.131) by taking the following gauge trans¬ 
formation: 


T = Miu 

As a result, we have 


( 2 . 1 ) 
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It is easily verihed that the gauge transformation fl2.1h does not change the mon- 
odromy matrix E^, k = 0,1. 

Let ('01,'02)^ be a vector solution of (12.21) . and set 


<!) = B 

then 0(A) solves the equation 
d?cj) 


dA2 


+ BC + A' - AB-^B' + -(B-'By - -B-'B" 

4' ' 2 


(2.5) 


( 2 , 6 ) 


where the derivatives are taken with respect to A, for instance. A' = 
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as X ^ +00, 


(2,7) 


Combining fll.lSp and fll.lbp with fll.Tp . we have 
v{x) = -| - \x"^F2{s) + ^ [F2(s) - 0] + O 
where -Fi(s) is defined in (II.8p . s = | — afeln and 
F2(s) = 06 ^“ - 


( 2 . 8 ) 


For X —)■ +CXO, substituting fl2.4p into fl2.6p . and in view of fll.7p and fl2.7p . a straightfor¬ 
ward calculation yields the following second-order equation; 


(P<p 

where 


= —X 


(2A-1)2 , gi(A,x) Q 2 (A,x) 


16A(A - 1) 
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X' 


3/2 


+ 0{x ^) > 0 :=-x^g(A,x)0, (2.9) 


Qi(A,x) = 


g2(A,x) = 


4A(A-1) _ 
1 

2A(A- 1) 


Aah + i — i { X 


H' 


FM + F^(s) - I A - i) F,(,)'^ 


2) H 
1 
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with 


as X ^ cxo. 


( 2 , 10 ) 

( 2 , 11 ) 

( 2 . 12 ) 


H{\,x) = ix'F (a - 0 + \Fi(s) - iF 2 (s) + O(x-F^) 

Here F\(s) and ^ 2 ( 5 ) are given in fll.8p and 02.81) . respectively. 

For large x, equation 02.9p has two turning points 

Xj = ^T x“^/^\/4a6 -h f (1 + 0 ( 1 )), j = l, 2 , (2.13) 

which coalesce with each other at 1 as x ^ cxo. The Stokes curves are defined as 

Re VA(A- 1) = 0. (2.14) 

In line with the idea of uniform asymptotics in [1], we define a number a such that 

■ ^ ‘ . (2.15) 


-TTia 

2 


= / {t^ — = / g^'^^(A, x)dA, 


' Ai 


and a new variable C, by 


(r^ — a!^)^/^(ir = f Q^F(j^^x)dr]. 

y. J X 2 


(2.16) 


Here and in 02.15p . the cut for the integrand on the left-hand side is the line segment 
joining —a and a. The path of integration is taken along the upper edge of the cut. 
With a and ( so chosen, the result in [U Theorem 1] applies. Thus we have the following 
theorem. 

































Theorem 2. Given a solution 0(A,x) of fl2.9p . there exist constants ci and C 2 such that 


/■2 2 _i 

( q(A x) ) ^ |[ci+o(l)]T)i.(e’"*/'‘\/^C) + [c2+o(l)]-D-j.-i(e""/‘^\/^C)| (2.17) 

as X ^ + 00 , uniformly for A on the Stokes curves defined in fl2.14l) . where D^{z) and 
D_^_i{z) are parabolic cylinder functions; cf. JF/l (12.2.5)], with 

V = —- + -ixa^. (2-18) 

2 2 ^ ^ 

3 The monodromy data for x — )► +(X) 


We proceed to calculate the two invariants Iq and Ji in fll.24l) - fll.25l] as a; —)■ +cxo. For 
our purpose, we need to clarify the relation between ( and A in fl2.16p . 

Lemma 1. For large x and X, 

^ixC - ^fxa^lnC = -|a + ^x + ^ + nab - ^ Inx + ^ ln(6e"^“) + o(l), (3.1) 


where 


= 


Aab + i 


X 


+ o\ - 

X 


as X ^ oo. 


Remark 2. Coupling fl2.18p and fl3.2p determines the approximate value 
u = 2iab — 1 + o(l) as x ^ cxo 

for the order of the parabolic cylinder function in fl2.17p . 

Proof. A straightforward integration on the left-hand side of fl2.16l) yields 


(3,2) 


(3,3) 


/ 2 2 
[T — a 


Y^'^dT= - {C(C^ - - o^ln (C +(c^ - +a^lna}. (3.4) 


Here, the cut for the integrand is again the line segment joining —a and a, and again we 
take the integration path along the upper edge of the cut. In view of fl3.4p and picking 
up the leading terms in fl2.16l) . for large C we have 

^C^-^a^ln(2C)-^a2 + ia^ln(o) + 0(o^C"^) =y' x)dr]. (3.5) 

To calculate the right-hand side of (13.51) . we split the integration interval, so that 


Q^/‘^{r],x)dr] = 


' A 2 


' A 2 


+ / \ Q^^‘^{'n,x)dr] ■.= I1 +I2 


(3.6) 
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where 


A* = ^ + (3.7) 

with T being a large parameter to be specified more precisely later. When A approaches 
1/2, it follows from fl2.1UI) and fl2.1ip that Qi ~ —{Aab+i) and Q 2 ~ —2(Fi(s) + F 2 (s)) = 
—4ae^*^, where s = | — a61n Fi and F 2 are given respectively in fll.81) and fl2.8p . Thus, 
in view of |Im{afe}| < A, we have 

= 0 ( 1 ) as a: —)■ cxo. (3.8) 

To approximate Ji, we make the change of variables 


replace Q\ by — (4a6 + i), and ignore Q 2 - Then for large x, we have 


h= Q^r],x)dr] = - _ 

J \2 ^ J \/4ab+i 

Aab + i Aab + i 

2x Ax 2x 


\/F — {Aab + i)^ (1 + o{l))dt 


ln(2T) 


4a6 + i 
Ax 


\n{Aab + i) + o [ — 
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(3.9) 


Here, as before, the cut for the second integral is the line segment joining —\/Aab + i and 
\/Aab + i, and the path of integration is taken along the upper edge of the cut. Taking 
T = —y/Aab + i in Ji gives (13.2p . 

When A is large, the three terms in the square brackets on the righthand side of 
fl2.1ip are of size o(x^^^) as x ^ 00 . Comparing this with fl2.10p . one can ignore Q 2 in 
h] cf. fl2.9p and fl3.6p . Accordingly we obtain 


-A r 


{2rj - 1 )^ 


(- 


Aab + i 


— i- 


277-1 


Jx* [ 1677(77 — 1) X \Ari{ri — 1) 877(77 — 1) 


if' 
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dr] 


. 277-1 
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00 

fAab + i i{2r] — l)H'\ 

A* 4^77(77-1) 

X (277 — 1)2 

V 4 8 h ) 


dr] 


(3,10) 


1 i 1 r" 

—A - — T — - - 

2 4 4 2x 


m{Aab + i) Aab + i 


, X 7 , 7X^/^ 

m —- H-In-^, 

T2 2 x ’ 


4x 4x 

with error term being o{x~^) + 0(A“^) + 0(T^x“^), where the derivative H' is taken 
with respect to 77 . Now setting T < substituting fl3.9p and fl3.10p into fl3.5p . and 
combining the latter with (13.2p . we obtain Lemma [H □ 


Now we have the relation between C, and A from fl2.16p for large A. We also need to 
establish the following relation for small A: 

Lemma 2. When A —)■ 0, for large x such that xy/X —)■ cxo and xCf^ —)■ oc, 

]^ixC - ]^ixa^ In C = ^a:\/A + + ^ ln(6e"2“) + o(l), (3.11) 

where given in fl3.2p . 
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Proof. Now, we assume that A —)■ 0 and x is large such that x\/~\ —)■ oo and —)■ oo. 

Let 


A^ 


1 

2 


— T^x 


1/2 


(3,12) 


where T* is a large positive constant to be specihed later. We split the integral on the 
righthand side of (I3.5p into three parts 


Q^^‘^{r],x)dr] = 


' X 2 


’ X 2 


' Xi 


+ / + / \ x)dr] Iq + Ii + I 2 . 


On account of 02.151) . we see that 

7ii 


* - 


Similar to 03.9p . we have 

rX* ^ rT 

h= Q^r],x)dr] = - _ 

J X 2 ^ J y/4ab+i 

Aab + i Aab + i 

2x Ax 2x 


— {Aab + iY (1 + o{l))dt 


Anh -\- ? 

ln(2r,) + -V->n(4oi + i)+t> 

4x 


(3,13) 


(3,14) 


(3.15) 


When A —f 0, for large a, from (l' 2 . 11 1 I we see that 

p Q2VVA-I),, 

A. iV^( 2 A - 1 ) 


(3.16) 


is of size o{'/Xx + o{x ^ ln(^)). Then, using the binomial expansion, the integral I 2 
in 03.13P is approximated as 


(2h-l)^ I 1 ( 


Aab + i 


2rj 


16r/(?7 — 1) X \Ari{ri — 1) 817 ( 7 / 

4a6 + i 


1 ) 


H' 

ll 


1 , 1 miAab + i 

-Ah-^ H—^- 

2 4 2x 2x 


IX 


dr] 

1/2 


Ax 


. X i . 

m —- H-In „. , 

Tf 2x ’ 


(3.17) 


with error term being o(x H- o(\/Ax H- o{x ^ln(^)) H- 0{T^x ^). Now choosing 
T* < such that o(a;“^ln(^)) = o{x~^), substituting 03.14p . 03.15p and 03.17p into 
03.13p . and combining the latter with 03.2p . we obtain 03.lip , thus completing the proof 
of Lemma [2l □ 


For large x, it follows from Theorem [2]that there are two uniform asymptotic solutions 
of equation 02 .9p (pi, and namely 


py{X,x) 


a 


Q(A,x) 



(3.18) 
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and 


= ( 1 ^) ’ . (3-«) 

the uniformity is with respect to A on the Stokes curves. 

We denote by tl/^ the {i,j) entry of \1/ and seek asymptotic solutions at inhnity and 
at the origin. By virtue of fl2.5p . we have 

^ + C20_._i), (3.20) 

= + (3.21) 

where Cj, j = 1, 2, 3,4 are constants to be determined by fll.181) . 

Similarly, we obtain 

+ C60-.-1), (3.22) 

+ C80-.-1), (3.23) 


where again Cj, j = 5, 6, 7, 8 are constants to be determined. 

For our purpose, we need to determine the asymptotic values of Cj, j = 1, • • • , 8. As 
a matter of fact, it follows from the hrst equality in fll.20p that 


= {Eo)u¥S + {Eo)2i¥S, (3.24) 

= (.Fo)i 2^S? + (^o)22^S?. (3.25) 


We obtain from fl3.20p . fl3.22l) - fl3.23l) and (13.241) that 


Cl — (F'o)iiC 5 + (F'o) 21 C 7 , C2 — {Eo)uCq + (F'o)2lC8, 


which in turn gives 

(A'o)ll _ _CiC8 — C2C7 

(770)21 CiCg — C2C5 

Similarly calculation leads to 

(770)22 _ _C3Cq — C4C5 

( 77 o)i 2 C3C8 — C4C7 

Eventually, substituting fl3.26p and fl3.27p into fll.24p yields 


(3.26) 


(3.27) 


7 o = 


(ciC8 — C2C7)(C3C6 — C4C5) 
(C3C8 — C4C7)(CiC6 — C2C5) 


(3.28) 


Now, we are in a position to calculate Cj for j = 1, • • • , 8. 

From [211 Sec.12.9], we have the asymptotic behavior of D^{z) for \z\ —?■ oo as follows: 


D^z) 


rsj 


z'^e 

z'^e- 





aigz e (-fTT, fir), 
arg^; e (^tt, fir). 


(3.29) 
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and the asymptotic behavior of D_y_i{iz) for \z\ —?■ oo is that 


D_y_i{iz) ~ 


VZn -■Ki{v+2M2 


on arg^; = — 
on argz = ^Tr. 


(3.30) 


For A on the Stokes line argA = —| and |A| —)■ oo, it immediately follows from fl3.ip 
that ~ if we take i = e^, then we have arg(C ~ 0. Therefore, arg(e”/^\/^C) ~ f 
and arg(e“^*/^\/^C) ~ — f for a; > 0. Taking —1 = e^\ from fl2.9p we have ~ 

\/2e-’r*/4 as |A| —?> oo. Since as |A| —)• oo, nsing the appropriate 

asymptotic formulas of D^{z) in fl3.29p . and in view of fl3.18l) and fl3.ip . we have 


01. ~ as |A| — oo. 


(3.31) 


Here use has been made of s = | — a61n |, and 


Similarly, by applying the appropriate asymptotic formulas of Dy{z) in fl3.29p . we obtain 
from fl3.19D and fl3.ll) that 

^-y-i ~ Boe~^^, as |A| —)■ oo, (3.32) 


where 

On the other hand, by virtue of fl2.ip . we obtain from fll.181) that 

(3-33) 

Thus, applying fl3.3ip and fl3.32p . a combination of fl3.20p . fl3.2ip . fl3.33p and fl3.34p 
implies 

C2 = C3 = 0. (3.35) 

Substituting fl3.35p into fl3.28p gives the simplihed version 

Jo = (3.36) 

CeO 

Our next task is to determine Cj for j = 5, 6, 7, 8. 

By using the result from [2l Prop.7.1], we have 

^^®Vj + o(1)) $ f -x\/A, e) (3.37) 

\nn + 0/ V2 / 
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as A —)■ 0 and x —)> oo, where Mi is defined in fl2.3p and 


$(x, Q) = X 


Je-i{x) Ji-e{x) 
Je+iix) J-{e+i)ix) 


(3.38) 


where J^ix) is the Bessel function of the first kind, and 0 ^ Z, as mentioned earlier. 

For argA = 27r, choosing A —?■ 0 such that x^/X —)■ oo and x(‘^ —)■ oo, by virtue of the 
asymptotic expansion of v{x) in fl2.7p and the behavior for Bessel function 


W = - 

TTZ 



COS ( z - -/ivr - -TT j ( /i 


-]sm(z--fi7i--7r]+0{z ') 


as 2 ; ^ cxo for I arg 2 ;| < vr; cf. 
fl3.37p as the following form 


, a straightforward calculation enables us to write 


^( 0 ) 


-Q f cos(|a:\/A - i07r - Ayr) - cos(|x\/A + i07r - Ayr) 

2\/F ^ I cos(|x-\/A — |07r — ^n) — cos(|a:\/A + |07r — ^ir) 


u 2 °"®, 


which gives us 


,(o) 




0 


4v^ 


+ i) ( 3 . 39 ) 


and 


^( 0 ) 

^12 


0 

40r 




(3.40) 


as x\/A — 00 . 

On the other hand, for argA = 27r, choosing A ^ 0 such that x\/~X —)■ cxo and 
t 00 , then we have from 03.111) that ~ \/A, which gives us arg^ ~ Ayr. Hence, 
arg(gW4y^(^) ^ ^ arg(e“’^*/'‘\/^C) ~ f for x > 0. From 02. 9 p we have ~ 

2A^/'^ as |A| —)■ 0. By using the appropriate asymptotic formulas of Dj,{z) in 03.29p . we 
obtain from 03.18P and 03.lip that 


~ 2A^/^x-^/^(C'oe-5^^ - Doei^^), 


(3.41) 


where 

Co = e-4-+i)/423-b-i/2^-i/2^ ^ (3.42) 

F(-z/) 

Similarly, by further applying the appropriate asymptotic formulas for D_^_i{—iz) 
in (I3.30p . we obtain from (I3.19p and (13.lip that 

~ ‘IX^I^x-^I^Fne^^^ (3.43) 
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with 


Fn = 


(3.44) 


Substituting fl3.4ip and fl3.43p into fl3.22p and using fl3.39p . we obtain 

0 


F1/22AV4^-1/4c5 


4:y/n 




(3.45) 


0 


~ + *)(e-5®"-3- + e^®’^+3^FoC'o"^). 

4v7r 


Similarly, substituting fl3.4ip and fl3.43p into fl3.23p and using fl3.40p . we obtain 

Ay/TT 


Therefore, taking ratios from fl3.45p and 03.471) . we have 

C7 

and a parallel result follows from 03.46P and 03.48P reads 


Cg 02®’^ 4^^ + e DqCq^ 


(3.46) 


(3.47) 


(3.48) 


(3.49) 


C 6 


= —u 


g 2®^ 4^^ + e2®^+4^FoC'Q ^' 


(3.50) 


where 


FnCn-' = 


— 1 _ ^y/irb ^TTg—7rab2~6ia6 


(3.51) 


“ “ r(l- 2 ia 6 )' 

Substituting 03.49P and O3.50p into 03.36p . we obtain the first invariant quantity for large 


X 


/n = 


ie"® - DqCq^ 

ie-TTiS _ DaCn^ ■ 


te — iyoO/Q 

For large x, in the same manner, the other invariant can also be calculated as 


(3.52) 


Ii = 


TTZ© _ 2^/7ra ^—FaboQiab 

r(l+2m6)^ ^ ^ 


te 


— 772© _ 2\/na 

r{l+2iab) 


^—'ivah2^Qiab 


(3.53) 


A combination of 01.26p and 03.52p with (13.5ip gives Ol.lOp . and a combination of 
fll.27p and 03.53P gives 01.9p . 
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